Isotropic Elastic Solid Ellipsoids 


39 


“ The Stresses and Strains in Isotropic Elastic Solid Ellipsoid 
in Equilibrium under Bodily Forces derivable from a 
Potential of the Second Degree." By C. Chree, M.A., 
Fellow of King’s College, Cambridge, Superintendent of 
Kew Observatory. Communicated by Professor W. G. 
Adams, F.R.S. Received March 2,—Read May 10, 1894. 
Abridged February 20, 1895. 

General Formulce. 

§ 1. Let the isotropic elastic solid ellipsoid, 

a~ 2 x 2 + b~ 2 y 2 + c~V = 1.(1), 

of uniform density />, be acted on by bodily forces whose components 
Vx, Q y, Rz are derivable from a potential 

y = i(P^+Q/+Rz 2 ).( 2 ). 

Let n denote the determinant 


3 h* + 26 V + 3c 4 , c 4 —y (b 2 c 2 -f c 2 a 2 + 3a 2 5 2 ), b i -—y(b 2 c 2 + 3 c 2 o 2 4* a 2 b 2 ) 

c 4 —y ( b 2 c 2 + c 2 o 2 + 3 a 2 b 2 ) , 3 c 4 + 2 c 2 a 2 + 3 a 4 , a 4 —y (3 b 2 c 2 + c 2 a 2 •+■ a 2 b 2 ) 

tf-y ( b 2 c 2 + 3c 2 a 2 -f a 2 b 2 ), a^y (Sb 2 c 2 + cW + a 2 b 2 ), 3a 4 + 2 a 2 b 2 + 3 Z> 4 

.. ( 3 ) 

and let its minors be n n , n m &c., where II 12 = n 2 i. 

Then with the notation of Todhunter and Pearson’s “History,” 
y denoting Poisson’s ratio, the stresses are given by 


-(iP,+8M + BT)J]> 

^ = & 2 [-(IQ/>+3N + L)^+aQ/'+N + L)(l-^ 

-(*<&>+N+sL)£]’ 

« = c 2 [-|(R P +L+3M) J-(iB*»+sL+M)j£ 


>■■■ (4); 


yz = 2L yz, 
YOL. LVIIT. 


+ (|R,+L + M)(i-Jj], 

^-V N 

zx = 2M zx, Xlj — 2 IS \xy 
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where L — [P a 2 {r] (b 2 + c?) n u 4- (yc 2 —a 2 ) IJ 12 + (^& 2 — a 2 ) ITi 3 } 

+ Q b 2 {( v c 2 -b 2 ) ILn + y (c a + a 2 ) n ia + (ya 2 -b 2 ) n 13 } 

+ Rc 2 {(^6 3 —c 2 ) ITn+^a 2 — c 2 ) II12 + 7 (n 2 + & 2 ) ni 3 }].. ( 5 ), 

while M and N are got from L by replacing the first suffix in the IPs 
by 2 and by 3 respectively. 

Denoting the displacements by a, /3, 7 , types of the 6 strains are 


«. = f; = * [(PV + M + N) [(iQp + N+L) ?> 2 

+ (iR / > + L+M)c 3 ] 

-'h(|P/'+M + N) a*- v [(iQp+8N + L) P 

hQRp + L + sM)^]} 

-£{QP P + M + 3N)a*- V [(iQp+N+L)b> >>... (6); 

+ (aE’/ 3 + 3 L+M) C 2 ]} 

-5{(-|P/° + 3M+R) a?- v [(iQ^ + K+3L) V 

+ (PV + L + M)c 2 ]}], 


— r]/3 dr-/ 4(1 + 7 ) 

R— L ^- 


I 


J 


where E is Young’s modulns, Poisson’s ratio. 
A type of the three displacements is 


« = J [(iPp + M + N) a 3 -7 QQ^ + N + L) & 2 

—7 (-|R/)4-L+M) c 2 

3 S {(Fp+M + N) a 2 -, (iQ^ + sN + L) & 2 

— f} (f E /9 + L + 3M) c 2 } 

-p{(iP/> + M + 3N) a 2 - 7 (|Q / , + R + L )& 2 

-7(AR / 9 + 3L+M)o 2 } 

{(iP^ + 3M+R> 2 -, (iQp+N + sL) W 

— V (^Rp + L + M) c 2 }] ... (7). 

The other strains and displacements can be written down from 
symmetry. 

For the elastic increment ca } in a principal semi-axis a , we have 
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oa'a = ~ [(iPp+M+N) a?- v {\ P (Qfe 2 +Rc 5 ) 

+ L(5 3 + c 2 )}] ... (8). 

§ 2. Let t x .... denote tlie component stresses across the tangent 
plane at x,y> z to a qnadric concentric with, and similar and similarly 
situa-ted to the bounding surface (1), for whose equation we take 


a“V + &"V + c"V = X... (9), 

then we easily find 

ts/ {|f(iP^+ M +N)« 2 } = h !{§! (iQ/»+N+L)&*} 

= *•/ {^ (4R/, + L + M)c2 } = 1_X .(W), 


where p* is the perpendicular from the centre on the tangent plane. 
Thus the resultant stresses across parallel tangent planes to the 
system (9) at the points of contact are all parallel, and their intensity 
varies as 1 —X. 

§ 3. If the ellipsoid be rotating with uniform angular velocity tv 
about the axis 2a, we have 

P = o, Q = R =2 tP ; 

while if it be gravitating, the force between unit masses at unit 
distance being taken as unity, 

P = 2.wpabc\ " 7 YrJ T" ' ' \TnS l_ — \7~T~a —TT.* • (H) # 

with symmetrical expressions for Q and It. If the ellipsoid be 
gravitating, and at the same time rotating about a principal axis, we 
have only to- add the respective values of P, Q.. It. Substituting the 
values of P, Q, It in the expressions for L, M, R, and inserting the 
consequent values of L, M, R in the formulae (4), (6), (7), we have 
the complete values of the stresses, strains, and displacements. 


Gravitating nearly Spherical Ellipsoid . 

§ 4. Denoting by /m the force between two unit masses at unit dis¬ 
tance, we may take 


P== 


;W r 


2 a 3 — V 2 —c 

5a 3 


with symmetrical expressions for Q and R. 
We thence find 


* Thomson and Tait’s Natural Philosophy ’ vol. 1 , Part IT, p. 47 . 

E 2* 
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2 7TjULp 2 a 2 


+${o+5y)0+zv)- 


^){( 7+5, )( 3-,? ) 


a 2 -6 5 


2 a 2 — b 2 —c 2 


(8—3 7} — 9^ 3 ) 


(2- v ) (3 + V ) | 


(2 + 3, + 15f)} 

a *' J 

+ S{ (7 + B,)(l + 3,)-^=^(a + 3, + 16,*) 

o i a 

(8-3,-SV)}] .(12), 


yz 


4 rrppSjZ 


16(1— v )0 + 5 v ) 


{ (7+£ 


’?)( 1 — 2 ''/) + 


2a 2 — b 2 —<y 


(1 —)?—4r/ 2 ) | 


(13). 


The other stresses may he written down from symmetry. 

With the notation of § 2, 

— 1—V 27rfip 2 x( S — y) f 2a 2 —b 2 ~ C 2 (2 — rj)(S + y) | 

* a. 15(1 — //) \ a 2 (7 4-5 97 )(3— y) J # 

It will he noticed that 


x = px' ( 1 —r 2 1r' 2 ), 

A- 

where r = v^O^-b^+s 2 ), V = v / (& ,2 +2/ ,2 4“ 2J,2 ) } 

a?', y\ z being the coordinates of the point where the radius vector r 
produced cuts the surface of the ellipsoid, and p being the perpendi¬ 
cular from the centre on the tangent plane at y', z'. "Near the 
surface we may put 

X — r 3 /r' 3 = 2(l — rjr r ), 


and so conclude that t x and the other stress components across the 
tangent planes to (q) vary approximately as the distance from the 
surface. 

As the stresses at the surface itself are of special interest in the 
event of any application to the earth, I shall briefly consider them 
for a spheroid in which b = a. The principal stresses are nn, tt, 00, 
directed respectively along the normal, the tangent to the meridian, 
and the perpendicular to the meridian. Using cylindrical coordinates, 
r f 0. z where 


r = \/(% 2 1ry 2 ), 0 = tan l y(x. 
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we get 
nn = 0, 

tt = 


1 


- 47 X c 




+ 


d 2 —c 2 1— v] — 4 ?? 2 | 


00 : 


15(1 — 7 ?) ~ 7/ ' a 2 7 + 5 rj /’ y 

—4 7riip 2 a? f a 2 — c 2 / 3r 2 \ 1 — t? - 4t? 2 1 

16 ( 1 =^) { 1_2 ’ + “^ _ ( 1_ tf) 7 + 5 7 / 


,(15) 


where p is the perpendicular from the, centre on the tangent plane. 

On the “ stress-difference ” theory of rupture an objection to the 
application to the earth of the results obtained by applying the 
elastic solid theory to a perfect sphere, is that the surface values of 

nn —00 and nn—tt would, for ordinary values of t?, be simply enor¬ 
mous .* This objection, however, ceases to hold when the earth is 
treated as incompressible and truly spherical, because tt and 00 then 
vanish, as well as nn. It is thus important to know what happens 
in the case of an incompressible material when the surface is slightly 
spheroidal. To do so, put y = ^ in (15), and we find 

nn = 0, 1 

^ | 

tt = 8t Tjip 2 (a 2 —c 2 )/285, y .. ( 16 )* 

00 = 8 7r/iyo 2 (u, 2 —c 2 )(l—3r 2 /a 2 )/2.85j 
Over the surface the maximum stress-difference, S, is the equa¬ 


torial value of tt —00, and is given by 

S = 87ryap 2 (n 2 — c 2 )/95 ...(17). 

Substituting for />, &, c values suitable to the case of a homogeneous 
“ earth,” we find that approximately 

S — 9‘4 tons weight per square inch.(18). 


This is large enough to show that even if the earth be supposed 
incompressible, the consequences of its mutual gravitation cannot 
safely be ignored. 

The strains and displacements in the general case of gravitation in 
a nearly spherical ellipsoid may easily be deduced from (6) and (7). 
From the expressions for the changes in the lengths of the semi-axes 
we get 


Sa 

a 


4 

15 


7 rpp 2 a 2 X 

“El 


1 — 2rj — - 


2 a 2 —h 2 — c 2 


2(l — 3 7 ? — 3?? 2 ) 

" 7+5^ 


} - ..(19), 


Sb Sc 
b c 


4 

15 


W 2 (^-c 2 ) 


(1 + t ? )(1 + 8 rj ) 

E (7 + 5-7) 


( 20 ). 


* See c Phil. Mag./ Sept., 1891, p. 247. 
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The principal axes all shorten in any material which is not very 
nearly incompressible. 

For absolutely incompressible material, 


ha _ 4 7r/LLpV 2a 2 — b 2 — c 2 

__ _ — — —- ^ " 


.-( 21 ); 


thus a principal axis shortens or lengthens according as the square 
on it is greater or less than the arithmetic mean of ttie squares on 
the principal axes. In a spheroid b — a, we find from (21), 

8a , be 4 7rfip\a 2 — c 2 ) 5 S /oox 

~7“ ir c“57~" E 6 S. ( 

where S is given by (17). 

Taking the numerical value (13) for 3, and for E the high value 
20 X 10 8 grammes weight per square centimetre, we should get from 
(22) for a spheroid the shape and size of the earth, a shortening of 
some 5 miles in equatorial diameters and a lengthening of some 10 
miles in the polar diameter relative to what these lengths would have 
been in the absence of gravitation. 


Rotating nearly Spherical Ellipsoid. 

§ 5. Suppose next that the nearly spherical ellipsoid rotates with 
uniform angular velocity w about 2 a. The values of xx, yy 9 xy , and 
yz are as follows :— 

t> 2\ r 

■j (7 + 5 i;)(3——5 if) 


- w"pa 


5 (1— y)(? ~i~5y) 2 


2 a 2 — b 2 — c 2 


2 (l -f"^) (6 — —-5 tj 




a 2 ~b 2 


b 2 


J 12 (7 + 5 V ) (3 - 6 V - 5 rf) - ~~ (39 + 16,,- 29'tf -10 ,f) 

““ (?* 1 
(9— Si] — 51 if —30 y 3 ) > 


,2 _/»2 


2 r a 2 —-b 2 

-< 2 (7 + 5t7)(3~67 ? — brf) - - (9- 8*} — 51 7j 2 — 30y 3 ) 


(39 + 16y-~297] z —10y 


'>}] 


(23), 





yy 
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(39 + 23 rj — 3 rf + 5 a; 3 ) - 


a -—C 


(6 + 13»; + 36i; 2 + 25 a; 3 ) 


} 

+ (l-|r) ( (7 + 5i?)(12 + V -5^) +-^— (13 + 9,, + + 5a; 3 ) 

a 2 ~c 2 , w 
-«r-(H-’?)(3+5»; 2 ) j 

- 2 (7 + 5ij) (4 + 7 rj + 5); 2 ) -+-j— (9 + 13); +27 a; 2 + 15-a; 3 ) 

o t ar 


- (6 — V — — 6^ 3 ) 1 


(24), 


— ix) £ pxy 


5(1— 77 ) (7-H 5 77) 2 

a 2 —b 2 

a? 

— ixPpyz 


(7 + 5 77 ) (3 — 6 ^ — 5 y 2 ) 


(27 +14^—9^/ a ) +~ 


- (3 +10 + 31 ry 2 -j~ 20 if) 


2 (7+5^)(4— 3 V ~5 V 2 ) 
2 a 2 -b 2 —-c 2 


yZ -o( 1-,;)(7 + 5 ,;) 2 L" 

+ ~~~ a a ( 3—4, ? —21 ’/ 2—10 ’? 3 )] 


(25), 


( 26 ). 


^ maybe obtained from i/y, and zx from xy by interchanging y with 
z and b with c. 

The redaction in length along the axis of rotation is given by 


(-&*/«) = 


2 tv 2 pa 2 


5E(7 + 5'^) L 


1 ^ 14 - 8 ^ 4 - 5 ^ 

2a 2 —b 2 —c 2 24 -f 223 t? 4- 282 y 2 + 95 if 
a 2 6 (7 + 


(27). 


It is thus greater or less than in a sphere of radius a , according 
as a 2 is less or greater than the arithmetic mean of the squares on 
the principal semi-axes. The tendency to shorten in the axis of 
rotation increases as the perpendicular diameters lengthen. 

For the increment in a perpendicular semi-axis 6, we get 


bbjb = 


r, 

5E(7 + 5tj) L 


8 —rj — 5a; 2 + 

+ 


6 2 + c 2 —2 a 3 25 + 16-a; —15 a ; 2 —10 a ; 3 


a 3 7 + 5 7] 

b 3 — C 2 (l + a;) (69+ 137 a; +70 a; 2 ) 


a‘ 


3(7 + 6i7) 


( 28 ). 
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From this we see that the “ extension per unit of length is 
greatest in the longer of the two principal diameters perpendicular to 
the axis of rotation* 


Very Flat Ellipsoid . 

§ 6. Supposing 2 a the short axis, we find, neglecting higher 
powers of cja and. b/a, 


xx — ip 


Pa 3 l-~ 


+Rc 2 


2(1-9) 


^ x y 


b 2 


0 ' 

-3 — 

G 


Pa 2 ( v b*-a*) +Qb 2 ( v a 2 ~-¥)-yRc 2 v (a*+b 2 ) 
4 a 4 -f 3 d 2 6 2 -r 4 F + yaPF 


zz — i pc 2 


P 1- 


Sx 2 


r 

6 2 " 


x 


a" (1-*— 


4f\ 

ar fc 3 l 

.... (29) 


a 3 6 3 


+ 3Jft(l- 


jr 

fe 3 ' 


0 Pa 3 (^£ 3 —a 3 ) + Q?> 3 ( V a 2 -h-) 
4 a 4 -{- 3 a 2 6 2 + 4 b 4 + >/a 2 6 2 


1 — 


2a? 2 2t/ 2 

l 2 *' 


^)]..(30), 


C = OTW { Pa 3 (yh*-a') ± Qfe 3 (yP-ft 3 ) + RcW+5 3 )} , 

J P V 4a 4 +3a 2 6 2 4-46 4 i~'i;a 2 6 2 ******** ^ * 


zx = •§■ pa?0 * 


-P + 


y Pc 2 
1 — a 2 


_ Pa 2 ( yb 2 - a 2 ) + Qb~(ya 2 — fr 2 ) + R c 2 y (a 2 + V) 
4 a 4 +3 a 2 6 2 -f- 4 6 4 + yd 2 F 


} ....(32). 


yy is got from xx, and yz from zx bj interchanging P with Q and 
a with b. 

The forces P, Q, P may occur independently, so the most important 
term in the coefficient of each has been retained. 

When P and Q alone exist, or the forces are perpendicular to the 

small dimension, zx and yz bear to xx, yy, and xy ratios of the order 
g : a, while zz bears a ratio of the order c 2 : a 2 . Thus for a first 

S ✓—S >*—S 

approximation the three stresses xx, yy , and xy alone need be 

retained. To a first approximation xx, yy, and xy are constant 
along any line parallel to the small dimension. These conclusions 


# Used here and subsequently in sense of total change in length of axis divided 
by original length, whether strain be uniform or not along the axis. 
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are so far in harmony with the assumptions usually made in theories 
of thin plates. It should be noticed, however, that dzxjdz is of the 

same order of magnitude as dxx/dx, and so of equal importance in 
the first body-stress equation. 

When R alone exists, or the forces are parallel to the small dimen¬ 
sion, zz at points in the interior is of the same order of magni¬ 
tude as xx, yy , and xy. 

If there be bodily forces both parallel and perpendicular to the 
small dimension, and the two sets be of the same order of magnitude, 
then to a first approximation the stresses due to the former set of 
forces may be neglected. 


Gravitating very Oblate Spheroid. 

§7. Taking /x as before for the gravitation constant, we have as a 
first approximation * 

P = Q = — 7 r^fipcja, R = — 4i7Tfip ........ (33). 

t As there is symmetry round the axis of z, we employ cylindrical 

N —S 

coordinates r, 0, 2 ?. The notation rr, 00, .... for the stresses 
explains itself. The displacements are u along r, and 7 parallel 
to z. The strains are 

s r = du\dr , Scf) — ujr, s z = dy/dz, "I 

, >■■■ (34). 

= 0 , — 0 , a zr = duIdzd^I dr J 

For the values of the stresses we find 


- 7r 2 p fac 


s +_v ( 1 _2 


<M> 


7 


11 + r A 


11+7 

— 2 7 rfip 2 c 2 (l—r 2 /a 2 — z 2 / c 2 ), 


( 35 ). 


ZT — 7 T 2 yp 2 


3 + 7 c 
Tl + v a 


r 0 = (j)z = 0 


The strains which differ from zero are 


# Thomson and Tait’s c Natural Philosophy,’ vol. 1, Part II, Art. 527. 





s r : 


3 up 2 ac (l — 7 f) f 

"i(ir+/l s+ ’- sl 


48 Mr. C. Chree. Isotropic Elastic Solid Ellipsoids under 

r s 1 1 




( 36 ). 




E( 1 X + /;) 


= 2 **rf C -zr ( liipSl+v} 

rr a E (ll+ 7 ) 

The displacements and dilatation are 
___ 7 r 2 pp 2 aor (l — rj} 


J 


E ( 11 -f y) 


( 3 + 7 — 


_ 2 7 r 2 jubp l ac (l — 2 7 ) 


h- CS7), 


E 


2 uc(l—2??) f r 2 l 

0 ^) {,+,- 2 ( 1 + ,)- t } 


Strictly, 00 , 0 r, and 7 , being of the order of terms omitted in 

—*N 

rr, s r , <fec., should be neglected ; but the record of their first approxi¬ 
mations may prove useful. 

rr, 00 and 00 are pressures at every interior point; is every¬ 
where a compression, and s z an extension, while s r is a compression 
inside and an extension outside the cylinder 

r = r 1 = a4HH^)/(l+^)} ... (38). 

The displacement u is everywhere towards the axis of the spheroid, 
and 7 away from the central plane, 0 = 0 . Every element is reduced 
in volume. 


Flat Ellipsoid Rotating about the Short Axis 2 c. 
§ 8 . Here 

P = Q = R = 0 . 

Putting for shortness 

4a 4 + 3u 2 & 2 -f4& 4 -b*/u 2 & 2 = IT ......... 


we find 


(39), 





Forces derivable from a Potential of the Second Degree . 49 

xx — a 2 pa % \_{a^-\-a 2 b 2 (X -+-?/) + 6 4 } (1— x 2 /a 2 ) 

~(i+*i)aY]+n', 

yy = iv 2 pb 2 [— (1 + 37 ;) b 2 x 2 

+ {a 4 + a 2 b 2 (1 +*,) + V}(l-y 2 /b 2 )-] -IT, 

^ >...(40). 

22 = — uy 2 pc 2 {a i +a 2 b 2 (l + ?/) + 5 4 } 

X {l -2x 2 la 2 —2y 2 lb 2 —z 2 lc 2 }-^U , y 
xy = — w 2 /j (a 4 ~2?/a 2 6 3 +& 4 ) asy-Ml', 

yz/yz = zx/zx = ^w 2 p{a i + a 2 b 2 + b i }--ll f J 

zz, it should be noticed, is negligible compared to yz and zx , and 
these in turn negligible compared to xx , yy and xy. The values of 

xx, yy, and xy bear to the corresponding first approximations in the 
case of an elliptic disc* of semi-axes a and b, rotating about its 
cylindrical axis with the same angular velocity tv, the common ratio 

3a 4 +2a 2 6 2 -f 36 4 : 4u 4 -f (3 + ?/) a 2 6 2 -f 4fr 4 .(41). 

This ratio also applies to the first approximations to A, s x , s y , s zi 
and a xy ; while the other strains are in both cases negligible compared 
to these. It applies further to Saja and hbjb. Thus the great 
majority of the results worked out for the rotating disc can be 
simply modified so as to apply to the flat ellipsoid by means of the 
following table which proceeds to three places of decimals. 


Table I.—Values of (3a 4 +2a 2 6 2 +3& 4 )~{4a 4 + (3 + ^) a 2 & 2 +4& 4 }. 


V 

bja 

0 . 

0*2 

0*4 

0-6 

0*8 

1*0 

0 


0*75 

0-748 

0-741 

0-734 

0*729 

0*72 

0*25 


0*75 

0-746 

0-735 

0-722 

0*714 

0*71 

0*5 


0*75 

0-744 

0-729 

0-711 

0*699 

0*696 


On either the “ greatest strain ” or the u maximum stress-difference ” 
theory of rapture, the limiting angular velocity in the fiat ellipsoid 

* ‘Phil. Mag./ July, 1892, pp. 70—100; see results (17), (18), (19), p. 75. 
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bears to that in the thin disc, of the same material and axes, a ratio 
which is the square root of the reciprocal of (41). 

If the limiting angular velocities be w for the ellipsoid, id for the 
disc, then the approximate values of u>ju f when y = *25, are as 
follows:— 

Table II. 

b/a = 0 0*2 0*4 0*6 0*8 1*0 

wjtd = 1*155 1*158 1*167 1*177 1*184 1*186 

In the case of the thin disc many very simple results were found 
for the radial strains; these apply without modification to the strain 
s r in the flat ellipsoid. 


Flat Ellipsoid rotating about one of its longer axes 2a. 

§ 9. Putting P — 0, Q = R = a) 2 , 
and employing II' as in (39), we find 

£= F fa'-v) ( i ■ -~F4) , 1 


4ar y 2 \ 1 

2x 2 2y 2 z 2 \ 

a 2 b 2 c 1 } 




X = ~k 0/a*—& s ) zx, 


xy = 0/»*—&*) xy 


>...(42). 


/-—s 

As b/a passes through the critical value the stresses xx, zx, 

and xy vanish and change sign; and at the critical value, to a first 
approximation, the only stress is 

yy = i"' V^'0 — x‘j a 1 — ifjb-) . (43) 


The stretches answering to (42) are 
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* E 


*/ + 


0 a?-rib 2 ) (I 2 ~- v a?) 


n r 


+ 3. 


(¥- v a 2 ) ( ¥ 


n' 


}(-!-!) ] 

)]• 




- _ ,"V'T r, »• !/- 

• ?/ “ 3 E Ll n' J \ U\ 


+ 


3(i 2 — 7)0?) I 6 2 2 a 5 


n' 




^ (44). 




x tf P b* v \- / „ , (^a 3 -& 2 )(a 2 +6 2 )U, * 2 / 

" 3 E Ll n' J \ a 2 &*, 


■ II' \ a 2 ^ ¥ J . 


For the increments in the semi-axes we have 


-f = -|^{2,+ (b 2 -,a 2 ) (2« 2 +,b 2 )/n'}r 

i^{2-(& 2 -,a 2 )(,a 2 + 2& 2 )/n'}, J..(45). 

* = {i-(& 2 -^ 2 ) (« 2 +^-)/n'} 

It is easy to prove 8b jb positive, and 8c-jc negative, for all values of 
7 ] and b/a; thus, of the axes perpendicular to the axis of rotation, the 
longer lengthens and the shorter shortens. The axis 2 a of rotation 
is always reduced. 

When the flat ellipsoid is spheroidal (— 8c jc) is greater or less than 
(— Saja), according as rj is greater or less than *25. 


Very Elongated Ellipsoid . 

§ 10. Treating ajc and b/c as both very small, we find 
xx . 4(1—jf) (3a 4 +2a 2 6 2 -|-3 ¥) 

= 2P^a 2 (l + v ) £{(1 — »7)(2a 4 +2 a?b 2 + 3 & 4 ) + i/a 2 & 2 } (l 

— {2 a 2 + 3 & 2 + »/(tt 2 —3& 2 )}?/ 2 J 

r 3//^ ¥ "1 

+2Q P VaHl+ v ){W+V)-V} j 
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- Bpo* [{ (1 ■- >f) + aW + V) - n (1 + 3 ,) a 2 b 2 } (l - ~j 

— aV (1+27]) (1—5)/) |- 2 

pin 

-{(l-f) (4a 4 + 3ffl 2 & 3 +46 4 )—'/(l + 3'7)a 2 5 2 }-^ • • ( 46 )> 

/y2 _i_ 7)2 

£2 (l — ?; 2 ) 4 (3a 4 42a 2 & 2 43 & 4 ) 4--— { (3 + t/) (3 a 4 4 2 a 2 & 2 4 3 & 4 ) 487;a 2 & 2 } 

L o - 

= Vpa^rf — 2{3a 4 42a 2 & 2 43& 4 — (a 2 — yb 2 ) (a 2 43& 2 ) —2 if a 2 (a 3 — 5 2 )} — 
~2{3a 4 + 2a 2 & 2 + 3 6 4 -(a 2 -7 / 6 2 ) (3a 2 4* 5 2 ) +2^ 2 & 2 (a 2 -& 2 )} -|j 
+ {3a 4 + 2a 2 5 2 + 3& 4 —-2 (a 2 —-^6 2 ) (a 2 4& 2 ) —nf(a 2 — & 2 ) 2 } ^1 — 

+ Q^& 2 ^ —2{3a 4 42a 2 5 2 43& 4 — (IZ—ya 2 ) (a 2 -\-3b 2 ) — 2r) 2 a 2 (a 2 — 6 2 )} 

—2{3a 4 + 2a 2 Z> 2 + 3& 4 — (b 2 — 7ja 2 ) (3a 2 +b 2 ) + 27j 2 b 2 (a 2 —b 2 )} 


-f{3a 4 + 2a 2 6 2 + 3J 4 -2(6 2 -7;a 2 ) (a 2 4& 2 ) - v 2 (a 2 -b 2 ) 2 } ( 1 


4JEVc 2 { 3 a 4 + 2 a 2 6 2 4- 3 4 v (5 a 4 4- 3 a 2 6 2 44 6 4 ) - fa 2 (3 a 2 + 6 2 ) ~ 5 f (a 4 - ?; 4 )} ^ 

- {3 a 4 4- 2 a 2 6 2 4 3 & 4 4 ?/ (4 a 4 4 3 a 2 6 2 4 5 6 4 ) - f& 2 (a 2 4- 3 ¥) 4- 5 v 3 (a 4 - 6 4 )} ^ 

c 

+ { (1 - ^) (3 a 4 + 2 a 2 5 3 + 3 & 4 ) + -4-Y 3 a4 + 2 « 2 & 2 + 3 6 ‘ + V (2 a 4 + 3 aV + 2 V) 


- v 1 (3a i + a 2 b 2 + 3b i ) -2 f(a i +b i )^jj .(47), 

xy . 4(1 — ^ 2 ) (3 a 4 4-2 a 2 & 2 4 -3 6 4 ) 

= /o#y[4(l4.^) P(^ 2 {7; (a 2 4& 2 ) —a 2 } 44(14?/) Q& 2 {^(a 2 4 & 2 ) — b 2 } 

4R{(1-j/ 2 ) (a 4 4& 4 ) 42 7 (l4 37;)a 2 & 2 }].. (48), 

^ f n 2 A~h 2 "1 

(l —t/ 3 ) 4 ( 3 a 4 42 a 2 b 2 43& 4 ) 4-£~~ {( 3 4t/) (3 a 4 4 2 a 2 6 2 4 3 fr 4 ) 487/a 2 6 2 } 

L c J 

= pzx 1 P~y{3a i +2a 2 b 2 i-3b i —4ib 2 (a 2 —yb 2 )—7} 2 (a 2 —b z ) ( 3a 2 -j-b 2 )} 

+ Q & 4/{3tt 4 +2a 3 6 3 +3& 4 -46 3 (6 3 - 1? a 3 ) -f(a 2 -¥) (3a 2 +& 2 ,i} 

_ R [ (l _ f )(3a 4 +2 a 3 6 3 +3 6 4 ) + ~^6 8 (3 a 4 +2 a 3 6 3 + 3 Z> 4 ) 

— r/ (3a 6 — 2 a 4 6 2 — a 3 & 4 — 2 //) — rfb" 1 (3a 4 +1 w'lr + 3 V) 


+ f (a 2 + i 2 ) (a 2 - 2 b 3 ) (3 a 3 + 6 3 ) 


.. (49). 
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yy and yz may be deduced from xx and zx respectively by inter¬ 
changing P with Q, x with y , and a with h. 

If P, Q, and It be of the same order of magnitude, the principal 

terms depending on them in the valnes of xx, yy, and xy are likewise 

of the same order; but in zz, yz, and zx, the principal terms in P and 
Q are only of the same order as the secondary terms in It. I have 
thus thought it best, in (47) and (49), to retain secondary terms in 
the coefficients of It, and to write the second approximation value of 
IT. If, however, R be zero, the terms on the left-hand sides of these 
equations with c 3 in the denominator should be dropped. 

When R alone exists, or the bodily forces are parallel to the long 

dimension, then, except near the ends of the long axis, zz is large 
compared to yz and zx, while these in their turn are, at most points, 

large compared to xx, yy, and xy. The hypothesis usually made in 
treating long rods, viz.:— 


xx z=. yy ■=. xy ~ 0 


(50), 


is thus approximately true, dxxjdx is, however, of the same order 

as dzxjdz, and dzxjdx of the same order as dzzjdz, so that the neglect 
of any of the nine differential coefficients appearing in the body-stress 
equations would be unjustifiable.* 

When the bodily forces are perpendicular to the long dimension, 

then—excluding special values of x, y, z — xx, yy, zz and xy are of the 

same order of magnitude, and are large compared to yz and zx. This 
result differs widely from (50). 

When P, Q, R are of the same order of magnitude, we may, for a 

/-s 

rough first approximation, neglect all the stresses but zz, and take 
zz = JR p(c z —z 2 ), 

i 

— «ar/ s 7 = — SyH = = 4 OV/E) (c 2 — 2 2 ), 

«jx = p/y — -i(R^/E) (c 2 — 2 3 ), 

> .(51). 

7 = !0fy/E>(e 2 —12 2 ), 

Baja = Bbjb — —JR/rayc 2 /!], 

Bcjc — ~Rpo 2 /E j 


The strain at any point is the same as in a long bar subjected to a 
tension at its ends equal per unit of section to the local value of zz. 

# See Todhunter and Pearson’s f History, 5 vol. 2, Part II, pp. 189—191. 
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To a closer degree of approximation we have, for the stress system, 


— zxjzx = —yzjyz = zzj{c 2 —z 2 ) = JR/> .(52), 

and for the strain system, in addition to s x , s y , and s Zi in (51), 


<rzx/zx — Oyvlyz = —JR pjn 


(53), 


where n is the rigidity. 

An analysis of (52) gives, on a plane perpendicular to z , a 
stress J Rp (c 2 — z 2 ) parallel to z , and a shearing stress —JR pzr 
along r, the perpendicular on the axis of 2 ; on a plane perpendicular 
to r a shearing stress, — JR pzr, parallel to z; on any plane containing 
the axis of z , no stress. 

Even in the general case with P, Q, R all existent, and secondary 
terms retained in the coefficient of R and in IT, we get for s & the 
simple formula 


s ff — 


P_ 

4E 


-V -X s - ( c ~ z ) 


+ R { (1 + (1 +7,,) (c*-**)-( 1 + 27 ) (rf+y 2 ) 

-i V (a 2 +6 3 )}]..(54). 


When the bodily forces are perpendicular to the long axis, the 
stretch parallel to that axis is thus appreciably constant over a cross 
section; these perpendicular forces tend to shorten or to lengthen 
the long axis according as they act outwards from it or towards it. 


Elongated Ellipsoid Rotating about the Long Axis 2 c. 

§ 11. Putting in (46) to (49) 

P = Q = w 2 , R = 0, 

we get 

— - w2/>a2 ( [{at+aW+V-r, (a i + b i )} (l -■ 


(1—7) (3a 4 + 2a 2 6 2 + 3& 4 )| 


- ti) prf 


( 1 — 7 ) (3a i +2aV + 3b i ) 


-( 1+2 7 ) ay 

^{a 4 +a 2 6 2 + 2 b *—7 (a 4 - 6 4 ) j * 2 
+ {2 a 4 + a 2 5 2 + 5 4 +7 (a 4 — 6 4 ) }y'‘ 


..(55), 


(a 2 +5 2 ) {(a 2 + 6 2 ) 2 - 7 (a 2 -5 2 ) 2 }(l- |)].... (56), 
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^ _ u 2 p{y(a 2 + b 2 ) 2 —a 4 —P}xy 

Xy ~ (1-t;) (3a 4 + 2a 2 6 2 +3& 4 ) V.. . . 1 h 


h?fiY)ZX 

4 ( 1 — rf) c 2 (3 a 4 -f- 2 a 2 b 2 + 3 5 4 ) 


[(a 2 + & 2 ) (3a 4 +2a 2 5 2 + 3& 4 ) 


-4b\a i -27 1 a 2 b 2 + b i )-f(a 2 -b 2 ) (a 2 + & 2 ) (3a 2 + & 2 )].. (58). 


yy and yz may be got from xx and zx respectively by interchanging 

x with y and a with b. To a first approximation yz and zx are 
negligible. 

Neglecting z 2 jc 2 in these formulae, we obtain results applicable to 
the central portion of the long ellipsoid; these results are identical 
with those I have previously obtained for a long elliptic cylinder,* 
the axes of whose elliptic section are 2a and 2b, rotating about the 
cylindrical axis. To deduce results for the elongated ellipsoid from 
those found for the infinite cylinder, we write 1—z 2 /c 2 for 1 in the 
constant terms in s x , s y , s z , and A, in the coefficient of x in a and in 
that of y in j3 ; and we multiply the expression for 7 by 1—^ 2 /c 2 . 

The strain and displacement parallel to the long axis are of special 
interest; they are 


s* = 


7 = 


4 


1 

4 


a>>? (V + & 2 ) /„ 

. 

E \ c 2 / 


,J PV (d+V) / , 

E l 1 3 c 2 /. 

..... (60) 


Denoting by 2 1 the length of the long elliptic cylinder, we have 


Sc/c = -i(o;V7/E)(a 2 + & 2 ) = KoZ/Z) ..(61). 


This enables the values of bc/c, for y ~ 0*25 and &/a = 0, 0*2, 0*4, 
0*6, 0'8, and 1, to be written down from Table XY, p. 159, of my 
paper on the elliptic cylinder. Table XYI of that paper gives values 
of caja and bbfb which apply unchanged to the long ellipsoid; while 
Tables XIII and XIY give its limiting angular velocity on the stress- 
difference and greatest strain theories. 

The formula 


Sa/a-Sb/h 


*V( 1 + 7 ) (a 2 -6 2 ){(l + 7 ) ( a 2 +b 2 ) 2 + a* + V} 

3 E (3 a 4 + 2 a 2 b 2 + 3 5 4 ) J 


shows that of the principal transverse axes the longer is that which, 
even proportionately, is most extended. 


* ‘ Phil. Mag.,’ Aug., 1892, formulae (80) to (83), p. 156. 
VOL. LYUI. 


F 
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Elongated Ellipsoid Rotating about a Short Axis 2a. 
§ 12. We have to substitute in equations (46) to (49) 

P = 0, Q — ft = w 2 . 


Putting for shortness 

4(1-V) (Sa i + 2a 2 P-\-s¥) = H" ... 

we get 

= -5rT{(W) (« i +3i')+(l-3 r 6,>¥}(l-| 


0J [J 0 ‘ 

yy ~ ~jpr 


'"■v nr pc 


— { (1 — f) ( 4 a ‘ * + 3 “ 2&3 + 6 V) — rj (3 + 5 rj) a?W) 

— {(1 — 7 3 ) (a 2 + 6 & 3 ) — 3'>](3 + 5i))a 2 }y' 

{ 5 a 4 + 3 a 2 6 3 -f- 3 b' 1 +3 Ija~b~ — rj 1 (5 a 4 — 2 a 2 b 2 + 3 6 1 ) } (.1 

-«*{(! + 3 r, + 4,rj) b 2 -\-2 (l — rj 2 ) a 2 } (j~ + 


( 63 ), 


, (64), 


j)]....(65), 


- { 3 a 4 + 2 a 3 6 2 + 3 6 4 + j; (5a 4 + 9 a 2 6 3 + 6 5 4 ) - </ 2 (3 a 4 - a 2 6 2 - 6 & 4 ) 


— rj 3 (a 2 —b 2 ) (5a 2 + 96 2 )} 

0 


— {3 u 4 4- 2 a 2 5 2 -f“ 3 5 4 +(10 u 4 + u 2 6 2 4-9 6 4 ) 4- >/ 2 (6 a 4 4* a 2 5 2 ~ 3 bS) 

+ 7 8 (a a -& a ) (5a 2 + 9^)}^//c 2 

4-{in"4-terms of order u 2 /e 3 and b 2 jc 2 } /i — JLA # . # .(66), 


yz = —±n) 2 pyz ... (67), 

ZX = — ...(68), 

= ~ir{(. I-'/ 3 ) (a 4 —3 ft 4 ) +2 7 (3 + 5 1? )a 2 6 2 } .(69). 


First approximations to the values of the stresses, strains, displace¬ 
ments, and increments of the semi-axes may be obtained by writing 
w z for ft in (51). 

To this degree of approximation it makes no difference which of 
the two short axes is that about which rotation takes place. The two 
short axes shorten to the same extent per unit of length, while the 
long axis lengthens. 

Py writing w 2 for Q and ft in (54), we get a very close approxima¬ 
tion to the strain s z parallel to the long axis.' 
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Application of Method of Mean Values. 

§ 13. Let t h ZI4a/3, t^ t 3 be the mean lengths of material lines 
parallel to the principal axes 2 a, 2fe, 2c, and let v be the volume of 
the ellipsoid. Then it is easily proved from, the results in my paper* 
on the mean values of strains, &c., that the elastic increments in the 


general case are given by 

i 0*/®) {Va 2 — 7 ) (Ofi *-f Re 2 ) }.(70), 

Zv/v = tVW*)(P^ + Q& 2 + Rc*) .(71), 

Mi/*i + %/*2+Ms/fe = &’/*.•. -..(72). 


In (71) A; denotes the “ bulk modulus.” 

For the gravitating nearly spherical ellipsoid, (70) gives 


«i/*i 


■A 71 '/'/' 2 ® 2 


1 ~ 2l > f , l-4^ 2a 2 -y— c 2 1 /-..-v 

E t “l-2« a 2 J ’ A 


Comparing (73) and (19) we see that in a gravitating perfect 
sphere 

< 5 £j/£i = r)a/a, 


or the reductions per unit length of a diameter and of the mean 
parallel chord are identical. 

For the increment in volume in the gravitating nearly spherical 
ellipsoid, (71) gives 

Sv/v = —47r/jp 2 (a 2 Fb 2 H-c 2 )/(45 7c). *.... (74).f 


For the gravitating very flat spheroid of § 7, we find 


Stifti — 7T 2 /ip%c(l 7j)l 5E, ) 

I 

dt 3 /t 3 = 27r 2 pp 2 ac't]l 5E, }> .... (75). 

Sv/v = —27r 2 /ip 2 ac(l —2?/)/5EJ 

For the general case of an ellipsoid rotating about 2a, we have 
btlih : 




£v/v = <o 2 /j (6 3 -j-c 2 )/(l5fc) 




( 76 ). 


# ‘ Camb. Phil, Soc. Trans.,’ vol. 15, pp. 313—337. 
f Agrees with formula (105), p. 335, ‘ Camb. Trans.,’ loc. cit. 


F 2 
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Isotropic Elastic Solid Ellipsoids. 

A material line parallel to an axis 2b, perpendicular to the axis of 
rotation, is exposed to two opposing actions. The components of 
“centrifugal ” force in its own direction tend to lengthen it, while those 
in the perpendicular direction 2 c tend to shorten it. On an average 
the former action will prevail so long as the mean dimension parallel 
to 2b bears to that parallel to 2c a ratio exceeding : 1 . 

Approximate Methods. 

§ 14. Suppose a body symmetrical with respect to the coordinate 
planes, of great length, 2 c, in the direction of z, to be acted on by the 
bodily forces whose components are Px, Q y, Rz. The stresses over 
the plane z = z must balance the force whose components are 

| JJ /3 Px dx dy dz , &c. 

Since every section perpendicular to z has its c.g. on that axis, it is 
clear the integrals vanish, which give the components parallel to x 
and y. If now we assume that when there is no rapid variation in 

the cross section, <r, as z alters, zz is large compared to the other 
stresses, ana is approximately uniform over <r, then we may take as a 
first approximation 



For an elongated ellipsoid this gives, as in (51), 

M = |R p(c z —z 2 ). 

Taking this as the sole stress, we obtain from the ordinary stress- 
strain relations the same values of s x , s y , s z as in (51), and thence by 
integration the correct first approximation values of the displace¬ 
ments. 

For a thin elliptic disc in which one of the axes, 2 c, of the elliptic 
section is very large compared to the other or to the thickness, we 
find as first approximations 

= ^R p (c 2 —S 2 ), 

«/* = Ply = J> . ..(77). 

7 = iRpz(c 2 -i2 2 )/E _ 

If, for instance, the disc be rotating about its thickness, i.e., the 
axis of the cylinder, then R = w z and 

ajx = Ply = —>V,0’-* 2 )/E, 

7 = 



( 78 ). 
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Micro-Metallography of Iron. 


As a third example, suppose the section perpendicular to the long 
axis uniform in shape and area, we then get 


& = |Rp(c 2 — z% 
u/x = p/y = -|R^(c 2 -/)/E, )> 


Rpz(c 2 -^)/E 


, n .(79). 


A special instance of this last case is presented by an elongated 
cylinder rotating about a perpendicular to its length through its 
centre.. Putting R, = to 2 , we have 


*/* = Ply — (« 2 — z 2 )/ e » 

7 = |wV(c 2 —I^O/E 


(80)* 


Comparing the several cases of rotation, we have an interesting 
illustration of how the effects of the “ centrifugal ” force increase as 
the.mean distance of the substance of the solid from the axis of rota¬ 
tion becomes larger. If the limiting angular velocity permissible in 
the elongated ellipsoid rotating about a short axis be taken as 100, 
then the limiting angular velocities in the thin elliptic disc and the 
elongated cylinder rotating about their short axes—the material, and 
the length of the long dimension being the same for all—are approxi¬ 
mately 87 and 71 respectively, both on the stress-difference and 
greatest strain theories. A caution must, however, be added that in 
bodies of such elongated form rotating about a short axis, a sudden 
change in the angular velocity may prove disastrous. 


“ Micro-Metallography of Iron. Part I.” By Thomas Andrews, 
F.R.S., M.Inst.O.E. Received December 15, 1894,—Read 
January 24, 1895. 

Secondary Micro-crystalline Structure in Metallic Iron . 

The term metallography appositely describes, in one word, that 
department of metallurgical science which deals with the accurate 
study and delineation of the ultimate formation or structure of 
metals, a knowledge of which is of the utmost importance. 

This development of the science of metallurgy is destined to prove 
of incalculable value in the study of the ultimate micro-structure 
of iron and steel. 

Experiments on the microscopic structure of iron and steel appear 

* This is more exact for a long beam of rectangular cross section than the result 
I obtained in the c Quarterly Journal’ for 1888, p. 29. 





